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Abstract 


An approximate method for calculating the noise generated by a 
turbulent flow within a semi-infinite duct of arbitrary cross section 
is developed. It is based on a previously derived high-frequency so- 
lution to Lilley’s equation, which describes the sound propagation in 
a transversely-sheared mean flow. The source term is simplified by 
assuming the turbulence to be axisymmetric about the mean flow 
direction. Numerical results are presented for the special case of a 
ring source in a circular duct with an axisymmetric mean flow. They 
show that the internally generated noise is suppressed at sufficiently 
large upstream angles in a hard walled duct, and that acoustic liners 
can significantly reduce the sound radiated in both the upstream and 
downstream regions, depending upon the source location and Mach 
number of the flow. 

1 Introduction 

Considerable effort has been invested in developing a new generation of su- 
personic transports. One of the primary requirements was that the aircraft 
be quiet enough to meet or even exceed existing noise regulations, and it was 
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decided that a mixer-ejector nozzle concept would be used to help accom- 
plish this objective. The idea was that a significant amount of the mixing 
noise would be generated internally within the nozzle, and could therefore 
be considerably reduced by using suitable acoustic liner designs. Data from 
recent tests using a prototype mixer-ejector show that the peak internal tur- 
bulence level is more than twice the external level. It is therefore important 
to develop prediction methods for this internally generated noise. A general 
theory based on Lighthill’s equation was developed by Goldstein and Rosen- 
baum [1]. Dill, Oyediran and Krejsa [2] extended this analysis to account for 
mean-flow refraction effects. However, both theories involve the solution of a 
complicated Weiner-Hopf problem, which can only be explicitly worked out 
for a slug ( or ‘top hat’) mean velocity profile. 

The experimental data suggest that the internal noise is of much higher 
frequency than the externally generated noise, in addition to being much 
more sensitive to nozzle geometry. Moreover, many of the most successful 
noise prediction schemes ( e.g. the MGB code) are based on high-frequency 
Lilley’s-equation solutions. The present study is therefore directed toward 
developing a high-frequency Lilley’s-equation solution that can be used to 
predict the internally generated noise. We suppose that the sound is gener- 
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ated by a superposition of statistically independent and acoustically compact, 
convecting-point quadruples and derive a formula for the high-frequency 
acoustic radiation generated by such sources when they are located within 
a semi-infinite, parallel-walled nozzle. In fact, we suppose that the mean 
flow is completely parallel, but allow the cross-sectional shape and veloc- 
ity profile to be arbitrary (as shown in Figure la), in order to account for 
nozzle-geometry effects. The only variation in the streamwise direction is 
due to the boundary condition change at the nozzle exit, which is allowed to 
have an arbitrary shape. Finally, an arbitrary (frequency dependent) acous- 
tic impedance boundary condition is imposed at the nozzle walls, in order to 
model an acoustically treated surface. The resulting solutions can then be 
superimposed to calculate the sound generated by an actual turbulent flow 
within a nozzle. The analysis can be used to guide the design of acoustic 
liners that may be required to absorb the noise radiated in specific directions 
or to design nozzle exit shapes that reduce the noise radiation below the 
flight path. 

Goldstein [3] developed a formula (equation (5.9) of that reference) for the 
high-frequency sound radiation from a convecting-point quadruple source in 
an arbitrary, transversely sheared mean flow. This result was later extended 
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by Durbin [4] to account for a general (not necessarily parallel) mean flow. 
These formulas involve a ray-spreading factor that multiplies the product 
of a source function-which describes the actual acoustic sources-with some 
Doppler factors that account for the local source and mean flow convection 
effects. The spreading factor accounts for the mean-flow variation along the 
path of the radiated sound and can be calculated from geometric acoustics 
or ray tracing. 

The present paper shows that Goldstein’s [3] formula still applies to the 
internally-generated noise and that only the ray-tracing analysis which is 
used to calculate the ray-spreading factor needs to be modified in order to 
account for the effect of the nozzle walls. This is demonstrated in Section 
2, where the notation is introduced and the Goldstein [3] and Durbin [4] 
analyses are reviewed in some detail. 

Three-dimensional ray tracing is fairly complex and somewhat difficult 
to implement numerically, but it was shown in Ref. [3] that the three- 
dimensional ray-tracing calculation could be reduced to a much simpler two- 
dimensional one for the doubly-infinite jet flow considered in that paper. The 
rays can then be found by solving a single second-order equation. In section 
3 we show that this can also be done in the present problem. The results 
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are applied to an actual turbulent flow in Section 4 and specialized to an 
axisymmetric mean flow in a round duct with circular exit plane in Section 
5, where some numerical results are also presented. Some conclusions and 
recommendations for further work are given in Section 6. 

2 Extension of Doubly Infinite Jet Solution 
to Account for Finite Nozzle Geometry 

For definiteness, we consider a unidirectional, transversely-sheared, parallel 
mean flow 

v = H7(x t ), p = p(x t ), c=c(x t ), p = constant, (1) 

with velocity v, density p, speed of sound c, and pressure p, exiting from a 
parallel-walled nozzle, as shown in Figure 1-a. The result which we obtain, 
however, is much more general and applies to more complicated flow configu- 
rations such as the one shown in Figure 1-b. Equation (1) is an exact solution 
of the inviscid, non-heat-conducting equations of motion for these configura- 
tions. x = {xj , x 2 , £3} denote Cartesian coordinates with Xi aligned with the 
direction of the mean flow, i denotes the unit vector in this direction, and 
x t = {x 2 , £3} denotes the transverse coordinate vector. The nozzle exit is 
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described by an arbitrary, three-dimensional curve C, as shown in the figure. 
The mean velocity U is assumed to go smoothly to zero at the generators of 
the nozzle wall and to remain zero beyond that surface. The analysis does 
not therefore account for forward flight effects, but can easily be extended to 
do so. 

Assuming that the ideal gas law applies, the linearized equation governing 
the acoustic propagation on this flow is [6] 

(2) 

where p now denotes the acoustic pressure fluctuation normalized by pc 2 , 


— = — II 9 
Dt~dt + U dxi ’ 


(3) 


denotes the convective derivative, and t denotes the time. T represents the 
acoustic source distribution and is given by 


r = 


7T Vf 

Dt 


2 v*y 


df 

dx^ 


(4) 


when this quantity is produced by a fluctuating force per unit volume. 
In the absence of temperature fluctuations, Lilley’s equation is obtained by 
replacing /, by the quadruple source distribution /,• = where u, denotes 

the velocity fluctuation within the flow. 
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Since the problem is linear, and the second term in (4) is negligible com- 
pared to the first in the high-frequency limit, the solution for an arbitrary 
force distribution fi can be obtained by superposition of solutions, say pc, 
to 



where u> is the frequency, x 4 denotes the source position, and S is the Dirac 
delta function. 

2.1 Review of Durbin’s high-frequency solution 

By using matched asymptotic expansions, Durbin [4] showed that the solu- 
tion to this problem is given by 

Pg = pg(x|x 4 ,lj) = (1 - Msx)$e tfc5 , (6) 

in the high-frequency limit 

k = u>/c 0 o — > oo, (7) 

where is the speed of sound in the region of zero mean flow, 

M = CZ/cco, (8) 
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S denotes the Eikenal, which satisfies the Eikenal equation 

(1 - M Sl ) 2 - (i) V = 0, (9) 

and 

s = {si,s 2 ,s 3 } = VS. (10) 

The solution to this first-order partial differential equation can be ob- 
tained by the method of characteristics by calculating S along the rays x(r ) 
, which are determined by the ordinary differential equations 


si = 0, 


( 11 ) 


Xi = Si 


1 


(I)' 


Uc 0 


c 2 


( 12 ) 


, i — 2, 3, 


( 13 ) 


Xi — 

« = j 

subject to the initial conditions at the source position x* that the initial ray 
velocity is proportional to the initial ray direction, say {cos //, sin n cos A, sin /i sin A}, 
i.e. that 


x* = 7 a{cos ju, sin n cos A, sin /i sin A}, 


( 14 ) 


where the dot denotes differentiation with respect to r, the subscript s de- 
notes quantities evaluated at the source position x a , and the proportionality 
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constant 7 4 is given by 


7s 


-2 



(15) 


As pointed out by Durbin [4], Eq. (15) without the subscript s, and with p 
replace by the polar angle 9 , relates the ray speed to the ray direction at all 
points along a ray. 

The amplitude function $ is given by 


$ = 


Pla s in P 


47TCCoo (1 - M,si) V PslJ 


(16) 


where J denotes the Jacobian determinant 


J = 


3 (xj , ^ 3 ) 

9(cr,/i, A) 


(17) 


with da = |<fx| denoting the distance along the ray. 

Once these equations are solved, the Eikenal can be found by integrating 
the equation 

S = s • x, (18) 


and the velocity fluctuation ug, corresponding to the acoustic pressure per- 
turbation pa, can be calculated from 

u G = c 2 s$e ifc5 /coo. (19) 
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The important thing to notice is that the derivation of these results is 
completely independent of any boundary conditions that are imposed on the 
surface S of the duct, and the termination curve C of the duct exit. The 
latter gives rise to the so-called defracted radiation which (Pierce [7]) is of 
higher order in frequency than the direct and reflected radiation and can 
therefore legitimately be neglected in the high-frequency limit-though it can 
certainly be important in the upstream direction. 

2.2 Modification of solution to account for the duct 
walls 

The conditions at the surface of the duct are accounted for by imposing 
boundary conditions on the solutions to the ray equations (11) to (13) at 
the point where the rays reach the boundary to produce a reflected wave, 
say {p+,u+}, corresponding to the incident wave, say {p_,u_}. (See, for 
example, [7].) 

The reflected wave is still given by equations (6) with (16), but multiplied 
by a constant reflection coefficient, say TZ. The Eikenal S is obtained by 
integrating (18) through the reflection. 
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Thus, the pressure and velocity on the boundary E are given by 

p G = $(l+Tl)e' kS , (20) 

and 

u G = ^(s_+7ls + )e* s (21) 

Coo 

(recall that U is assumed to be zero at E). 

The usual impedance boundary condition for a locally-reacting surface 
involves only the normal component of the velocity u G , and therefore only 
the normal component of the propagation vector s. This condition is usually 
expressed in terms of an impedance, say Z (which can, in general, be a 
function of the frequency u;,), as 

Z = pc?-M^ for x on E, (22) 

u G • n 

where n denotes the unit normal to E. Moreover, the normal component of 
s changes sign, i.e. 

s_ • n = — s + • n for x on E. (23) 

Substituting this along with (21) into (22) yields the following expression for 
the reflection coefficient 

n = ( 24 ) 

V - 1 
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where 


77 = - (s_ • n) C, (25) 

where £ = ZfpooCoo is a normalized impedance. Notice that Z — » 00 and 
7£ — » 1 for a hard wall. 

On the other hand, the tangential component of s, namely Si, remains 
unchanged by the reflection, and it therefore follows from (11) that 

Sj = constant, (26) 

which is equal to the far-field value of this quantity for any ray that propa- 
gates to infinity (which are the only ones we are interested in here). In this 
region (where the mean flow is zero), the acoustic rays are straight lines and 
are therefore given by 

x = x* -f /^(cosflocsinfloo cos^o^sinfloosin^oo) , (27) 

where R can be taken as the distance between the source point and the ob- 
servation point, and and 4 } <x> denote the far-field polar and circumferential 
angles, respectively, shown in Figure 1. 

It therefore follows from Eqs. (11) to (13) and the Eikenal equation (18) 
that R = 1, and that 

Si = cos^oo. (28) 
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The Jacobian determinant (17) becomes 


J = R 2 sin 0 o 


B{0 oo,4>oo) 

d(fx,X) ’ 


and it now follows from (12), (14), and (28) that 


7 a cos fJ, 




n , U S C oo 

cos Voo H zr— 


/ oo I _ 0 

c: 


Eliminating 7 * between this and equation (15) shows that #oo depends 
only on /i, and not A, and that 


3 / c 3 \ ^ sin 0oo 
s Vcoo/ sin d\i 


Inserting this into (16) and using (29) shows that 


(1 - M 5 cos0oo)^ Id^oo)’ 

where we have used the ideal gas law to obtain this result and put 


n = f[Ki, 


where the denote the individual reflection coefficients for each of the m 
reflections that the ray undergoes before leaving the duct. Also, it follows 
from (13), (18) and (28) that 


S = (x x - Xi“) cos 0oo + 5o(xt|x t 4 ). 
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2.3 Application to moving point source 


As indicated in the Introduction, the sound radiated by an actual turbulent 
flow can be calculated in terms of the pressure field p generated by a su- 
perposition of point quadrupole sources moving downstream with the mean 
flow. We therefore consider the source distribution 

r = (35) 

where U c denotes the convection speed of the source whose strength is Qij- 
The corresponding acoustic field can be calculated from the fixed source 
solution pg by superposing Fourier components and using the Green’s for- 
mula ([3]) 

p = t 1 / £> o(x|y ’ - x ‘ ” lu ^ iT ^ 

(36) 


Integrating by parts to transfer the derivatives from the source term to the 
Green’s function, and carrying out the integrations with respect to y ( and r 
gives 


p = 


Qij - 


iw,t 


U* 




d 2 


— Pg(x|x s , u)dx i a du. (37) 


2tt U~ J-coJ ' dxi’dxj 1 

For clarity, we begin with the case where only a single ray reaches the 
observer. The result will then be corrected for multiple ray effects in a 


15 



relatively obvious manner. Inserting equations (6), (28) and (34) into (37), 
and using the fact that (at lowest approximation) the partial derivatives 
operate only on the frequency-dependent terms in the exponent, we obtain 

p = (1 - M cos^ 00 )e- iu " 1 ^- f f 

x a fl ■ ^ e^ t(l| - Zl,)cosg °° +5ol dx 1 J da;. (38) 

oxi dxj 

Then, carrying out the integration, first with respect to xi* (to obtain a 
J-function), and then with respect to w, shows that 


(1 A/ COS 0 oo )<7 x <7 3 Q i j 2 ( - _ ” «■ ■ o , ^ y (n cos Oco +5 0 - g,, t) 

(1 — M c COs0oo) 


where we have put M c = U c /c <*,, k , = uja/c^ and 


cos Oca 

1 — M c COS #oo 


1 — M c cos 0oo <9x,' 


-, for 2 = 2, 3. 


Then it follows from equation (32) that 

| i2 k'lQwl it j dA 

^ (47rf?) 2 c^(1 — M s cos flco) 2 (1 — M c cos ^) 2 d</»oo 

in the far field where M — 0, which, except for some minor notational 
changes, and the inclusion of the reflection coefficient it, is the same as 
equation (5.9) given in Ref. [3]. The normalized wall impedance, £, which 
appears in this equation through the reflection coefficient, it , must be eval- 
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uated at the actual (or observation) frequency 


U) = U>,j (1 — M c COS^oo ) , 


(42) 


and not the source frequency u a . 

It is convenient to allow the transverse orientation of the quadrupoles to 
vary with source position. This amounts to changing the orientation of the 
x 3 coordinate system or, equivalently, referencing the angle A to a different 
angle, say A 0 (x 3 ). Then it follows from the results given in Ref. [3] that cr 2 
and cr 3 are given explicitly by 


cr 2 = 


_ — g 5 cos (A - A 0 ) 


1 — Mr cos 9 o 


03 = 


—q, sin (A - A 0 ) 

1 — Mr COS 9oo 


(43) 


where 


9 = 


\ 


(1 — M cosOooY 


— COS 2 One- 


(44) 


{cl Coof 

When multiple rays (which we individuate by a superscript in parenthe- 
ses) reach the observer, the far-field pressure is given by the somewhat more 
complicated formula 


Irf 


b 3 QtjQhDm 


(4nR) 2 (1 — M s cos ^oo) (1 — M c cos ^oo) 


2 y 


(45) 


where ★ denotes the complex conjugate, and the dependence on the transverse 
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source coordinates r 3 ,<f> s , and the emission angle X enters through 

i ^(S° ( " ) -So<"' ) ) 


<9A(") 

d>S m ) 

d<t>oo 

d(f > oo 


D ijkl = Y1 

n,m = l ' 


where k denotes the number of rays reaching the observer. 


(46) 


3 Reduction of Order of Ray Equations 

Goldstein [3] introduced the two-dimensional ray distance S defined by (see 
equation (2.19) of that reference), 

|dx t I 


dS 


= 1 . 


(47) 


It follows from equations (14), (15) (generalized to arbitrary position 
along a ray), (44) and (28) that r is related to S by 

dS 


dr 


= <?• 


(48) 


Equations (13) can then be combined to obtain the second-order system 


d dx t 

dS Q dS ~ tQ ’ 


(49) 


where V f denotes the cross stream divergence. This is the same as equation 
(2.23) of ref. [3], where it is shown, by introducing the polar coordinates 

4> = tan -1 (ar 3 /x 2 ) , r - y/xf + x 3 2 , (50) 
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that it can be reduced to the single second-order equation 

1 d r 2 q dq 

ld4>~T = 50’ 

where 

H r2+ ® ■ 

which is to be solved subject to the initial conditions 


r = r a ; — = r 3 cot (A - cf> s ) , at <j> = 4> s . (53) 

acp 

Inserting equations (13) and (48) into the boundary condition (23), using 
(50) and taking (j> as the independent variable, shows that the appropriate 
boundary condition for equation (51) is 


-r + tan(0-/?) = - r - tan (4> - (5) , for x on S, 

1+ d(p 1— d<p 


where we have put 


n = {cos/?,sin/?}. 


This boundary condition must be imposed on all rays reaching the cylindrical 
surface containing the duct wall whenever 


xi ^ x c [x2i *^ 3 )) 
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where ii = x e (£ 2 *X 3 ) is the equation for the termination curve C. Xi can be 
calculated as a function of <f> along the ray by inserting (28), (48), (50), and 
(52) into (13) to obtain 


dx i 
d(j) 



COS 0oo + M 


(57) 


4 Application to Sound Radiated by Actual 
Turbulent Flows 


Equation (45) can be used to calculate the power spectral density of a spec- 
tral distribution of sources of band width At * ) 3 by putting QijQli equal to 
. However, pressure spectra are measured per unit observation 
frequency Au; 

Au>« 

Aw = - A ;-— r , 

1 — M c COS Voo 

(see Eq. (42)) and it therefore follows that the directivity of the spectra at 
constant source frequency u 9 (due to a source at x*) is given by 



— I pi 2 
Aw m 


ks 4} $ijkiDijki 


(4nR) 2 (1 - M 4 COS0oo) 2 (1 - M c C OS0oo) 


d A 


dcj > o 


(59) 


This result can now be used to calculate the sound emitted by an actual 
turbulent flow by assuming that the turbulent eddies behave like compact 
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sound sources, and using Lilley’s equation to show that the spectral source 
strength 'J'.-jw is related to the fourth-order, two-point, time-delayed correla- 
tion function of the turbulence 


R ljkl {x a ,Z,T) = - < U 'j 


(60) 


in the usual way by 

Vijki = J J e~ iW3T Rijki{x s ,Z,T)dSdT, (61) 

where the single prime indicates that the quantity is evaluated at the po- 
sition and time (x J, ,t), the double prime indicates the position and time 
(x s ",f + r), 

£ = x s " - x s/ - i U c t, (62) 

and 

x-={*r,s«+*r) ,s (<+<’)} ( 63 > 

denotes the mean position of the source. 

Since the sound field is always produced by a distribution of sources rather 
than by a single point source, the final result will involve an integral of Eq. 
(59) (and consequently of Eq. (46)) over the transverse source coordinates r s 
and 4> s . Then, since S^ n) is a function of these coordinates, the contribution 
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from the cross-coupling terms in (46) will be smaller than the contribution of 
the m = n terms by a factor of (at least) which in a strict asymptotic 

sense is negligible in the high-frequency limit. However, the zero-mean-flow 
computations of Boyd, Kempton and Morfey [8] suggest that the asymptotic 
convergence may be relatively slow for sources close to the wall (which result 
in small values of S 'q — Sq in Eqn. (46)), and that the interference effects 
may not be insignificant even at relatively high frequencies - particularly 
at small angles to the downstream axis where the sound field in expected 
to be maximal. However, the turbulent flows, which are of interest here, 
will probably introduce significant random fluctuations in the phases of the 
disturbances, which will tend to uncorrelate the pressure fluctuations corre- 
sponding to different ray paths. We therefore feel that it is best to neglect 
the interference effects, which amounts to replacing Eq. (46) with 

D ijU = E <r>V n V n W B) |£ (n, | 2 • (64) 

Since the fourth-order correlation tensor is very difficult to measure ex- 
perimentally, or even calculate numerically, it is usual to assume that the 
turbulence is quasi-normal, and, consequently, that Rijki can be expressed as 
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the product of second-order correlations ([9], [10]) 


Rijki = RxkRn + RuRjk. ( 65 ) 

In order to simplify this further, Goldstein and Rosenbaum [11], Kerschen 
[12], and, more recently, Bechara et al. [13] and Khavaran [14] , assumed that 
the turbulence is axisymmetric about the direction of the mean flow. The 
analysis given in Ref. [1] (see also [2]) then shows that 

wnc = | of/ J €~‘ w,T QoiidT — 2(*T2<T3) 2 J /e" 1 ”’ 1 (Q, 3 - Qn) d(dr 

+ j J e-^Q^dr, (66) 

where we have dropped the superscript (n) on the <7,-, and 


On 

O 12 = 0 13 
Ql% = Q 33 
O 23 
Oo 


= Rh - Rn 

= Ri 2 + R\iRi2 
= R22 R12 

= R^2 — R-23 

= Rh 


are symmetric in their indices. 


(67) 
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5 Application to Round Duct with Axisym- 


metric Mean Flow 


Goldstein [3] showed that Eqs. (51) - (53) can be solved analytically when 
the mean flow is axisymmetric. A similar procedure can be used to obtain 
an analytical solution to the present problem, but it is probably easier to 
solve it numerically. However, it is important to notice that, in this case, the 
resulting solution, whether obtained analytically or numerically, will depend 
on A, 4> and <f>s (where (f> 3 is the circumferential angle of the source point) 
only in the combinations 4> — (f> s and A — since the coefficient q in Eq. (51) 
is independent of <f > , i.e., <f> appears only as an independent variable. This, in 
particular, implies that A — 4>, is a function of (f > and 9^. Moreover, 
calculations of the ray trajectories for sources located within the nozzle show 
that A is a discontinuous, multi-valued function of 4 > due to the sudden 
change of boundary conditions at the nozzle lip. This is illustrated in figure 
2, which is a plot of ^ vs. A for the indicated source location. Thus, even 
though 0oo is necessarily a single-valued function of A, the figure shows that 
the converse is certainly not true. It follows that 


p fj sin p (A< n > - <j> s ) cos 9 (A< n > - <f> t ) |ft< n) | 2 


dA< n > 

dftoo 


d<t> 3 
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d<f> s (68) 


= y "'f sin p (A*") — <£*) cos’ (A*") — <£,) |7^ n ) 

n=l JO 


d( A< n > -</>,) 

d<f>, 


-P 

Jo 


sin p A cos 9 A 




|2 - 


d\ for p,q = 0, 1, 2, 


since, for a given r 3 and #<*>, It depends on </> — <f>s only through A — <f> s . 
The sum in Eq. (68) must be taken over all Al"l values corresponding to any 
given value of <f ) <*,, in order to account for all of the rays reaching a given 
observation point. 

If we now choose the reference angle A 0 in Eq. (43) to be equal to 0 S , 
the quadrupole sources will have the same orientation relative to the radial 
direction for all <£ a , i.e. the quadrupole source distribution in Eq. (59) will be 
axisymmetric when Vfhyjy is independent of (f> a . Then, since equations (43), 
(44), (64), and (68) show that the entire <f>oo dependence in (59) is of the form 
(68), it follows that the sound field |p| emitted by a ring of uncorrelated, equi- 
strength quadrupole sources with radius r a , and the same orientation relative 
to the radial direction, is independent of the circumferential observation angle 
<£oo, i.e. it is axisymmetric. 

When |^.| = 1 (i.e. for a hard-walled duct), it follows from equations (43), 
(44), (66) and (68) that 


£ /%< n V n W n W n) *i;*, |^ (n f 

n= 1 J ° 


dA< n > 

d<t> oo 


d(f> s = 
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( 69 ) 


(i - [ v / / ■ e ' i “ ,T ( 8<5 ° + + ^ 

+2 cos 2 O^q? j J e 1U,,T + Qo) d£dr + 

+ cos 4 ^oo J J e IW,T (On + Oo) d£dr , 

For isotropic turbulence Q,j = 7Q 0 , and it follows from Eqns. (43) and 
(44) that Eq. (69) is independent of 6^ when M 3 = M c . This means that 
the sound radiated by a ring source in a hard-walled duct is not only inde- 
pendent of the mean velocity profile within the jet, but is also unaffected by 
the presence of the duct when all rays reach the far field. Of course, this 
result only applies when the phase cancelation between multiple rays can be 
neglected. Also, since r s d<f> a is the element of arc length, the total sound 
radiated by the ring source will be directly proportional to the radius r„. 


5.1 Numerical results 


Results for the directivity patterns due to a ring source within a round 
duct were computed for a constant mean speed of sound, c = Coo, and mean 
Mach number profiles of the form 


M(r) = Mo 



(70) 
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where Mq is the centerline Mach number and the parameters a and 6 are 
used to control the profile shape. 

The source terms in Eq. (66) were evaluated using the relations given by 
Khavaran [14] for axisymmetric turbulence. The anisotropy is characterized 
by the two parameters u\lu\ and where and u\ are the streamwise 

and transverse mean square turbulent velocities, respectively, and L\ and 
L 2 are the corresponding correlation lengths (see ref. [14]). Values for the 
anisotropy parameters of = 0.6 and Z^/Ti = 0.5 were used in the 

calculations. 

Figure 3 shows the results for the far-field directivity 

2 |p | 2 (47 rfl ) 2 g, f71 x 

Au, k^nn ’ 

vs. far-field polar angle, 0^,, for a ring source at r„ = 0.75r o ,Xi a = —2.0 r 0 , 
where ro is the duct radius, and a centerline Mach number of 1.5 with 
a — 0.1,6 = 6, for a hard- walled duct and a soft-walled duct of various 
impedances. 

For sufficiently small far-field polar angles outside the zone of silence, all 
rays emanating from the source reach the far field and, for the perfectly- 
reflecting, hard-walled duct considered here, the duct has no effect on the 
far-field sound. At far-field positions beginning in the upstream quadrant 
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( i.e. 0co > 7r/2 ), however, some of the rays become trapped within the 
duct, causing the sound pressure levels to be reduced at these angles. The 
hard-walled duct, therefore, only effects the sound field at sufficiently large 
angles to the downstream axis which, in fact, lie in the upstream quadrant 
as indicated in the figure. Since the number of rays reaching the far field 
rapidly decreases as 6^ — > it, there is a sharp drop in the far-field sound. 

However, the soft-walled duct starts to effect the sound field as soon 
as wall reflections begin. Since an increasing number of rays reflect (an 
increasing number of times) off the walls as the polar angle increases there 
is a substantial decrease in the far-field sound relative to the hard-wall case. 
The wall impedances £ = (1,-1) and ( = (2,-1) axe seen to reduce the 
peak noise level by nearly 5 dB, relative to the hard wall case. The results 
suggest that the magnitude and phase of the normalized wall impedance can 
significantly effect the peak sound level, and a detailed parameter study to 
find the optimal value should be carried out. 

Figures 4 and 5 show the effect of the source position on the far-field 
sound. The rays undergo fewer wall reflections when the source is closer to 
the nozzle exit (Fig. 4), and the acoustic liner therefore provides less noise 
suppression. When the source is closer to the duct centerline (Fig. 5), all 
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rays exit the duct without reflecting off the wall when the far-field polar 
angle is sufficiently small, and the acoustic liner has no effect on the sound 
field. Wall reflections start to occur when the polar angle is increased, and 
the liner reduces the far-field sound, but only by a relatively small amount - 
again due to fewer wall reflections. 

Figure 6 illustrates the effect of centerline Mach number on the liner 
effectiveness. At the subsonic Mach number (A/(0) = 0.9) for which this 
result was obtained, a wall impedance of ( = (1,-1) again reduces the peak 
sound pressure level by about 5dB, but produces a much larger reduction 
than the previous (supersonic) case at large upstream angles. 

6 Conclusions and future work 

It was shown that the high-frequency Lilley’s-equation solution devel- 
oped in Ref. [3] for a doubly-infinite, transversely-sheared mean flow also 
applies to the noise generated internally within a nozzle, provided appropri- 
ate boundary conditions are imposed on the ray trajectories at the surface 
of the duct and a suitable wall impedance factor is included. 

By assuming the turbulence to be axisymmetric about the mean flow di- 
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rection, a simplified expression for the far-field sound radiated by a turbulent 
flow within the nozzle was derived. 

The analysis was applied to the case of a round duct with an axisym- 
metric mean flow, and it was shown that a hard-walled duct has no effect 
on the far-field sound radiated at polar angles sufficiently close the duct axis 
(but outside the zone of silence). The numerical results show that the duct 
cuts off some of the rays for polar angles in the upstream quadrant, and that 
acoustic liners can significantly reduce the far-field sound but their effective- 
ness depends upon the wall impedance, source position and mean flow field. 
The analysis can be used to carry out detailed parametric studies to find the 
optimal wall impedance, acoustic source distributions, mean profile shape 
and nozzle geometry for a given application. 

The ray acoustics solution has the advantage of being applicable to nozzles 
of any shape and any mean velocity profile (see Fig. 1). The high-speed 
civil transport was expected to use a rectangular mixer-ejector nozzle with 
a very complex mean velocity profile and acoustically treated walls. Future 
work will evaluate the ray acoustics solution for this geometry and make 
comparisons with some recent test data. 

This paper only considers the sound produced by the first source term in 
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Eq. (4) - the so-called ‘self-noise’ term. While this term is asymptotically 
large compared with the second (or ‘shear-noise’) term in Eq. (4) in the 
high- frequency limit, it may be necessary to include the latter in order to 
obtain agreement with experimental results, particularly in the downstream 
quadrant (see, for example, Khavaran [14]). 

This paper also does not address the diffracted radiation produced by 
acoustic rays striking the duct lip. It too is asymptotically small compared 
with the direct and reflected sound (Pierce [7]), but can still be of significance 
at the upstream polar angles, where most of the direct or reflected sound is 
cut off by the duct. This is currently being investigated by Wundrow and 
Goldstein [15], who plan to develop a computational algorithm incorporating 
the diffraction effects into the present analysis. 


The authors would like to thank Dr. James Bridges of NASA Glenn 
Research Center for providing values of the wall impedance and Dr. Abbas 
Khavaran of Dynacs Engineering Co. for information on the mixer-ejector 
nozzle configuration and flow field. 
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Inner duct 










Figure 4: Far-field spectral directivity for Mach number profile (70) with a = 
0.1,6= 6,Af(0) = 1.5, and source position r t = 0.75ro,ri' = — 0.5ro for hard- 
wall duct (solid) and soft- wall duct with £ = (1,-1) (dashed). 






List of Figure Captions 


Figure - 1 (a)Flow configuration, (b) Example of more complex con- 
figuration to which analysis applies. 

Figure - 2 Far-field cirumferential angle vs. initial cirumferential angle 
for Mach number profile (70) with a = 0.1,6 = 6, M( 0) = 0.9 and 
source position r s = 0.75r o ,Xi a = -0.5r o , and 9^ = 3tt/8. 

Figure - 3 Far-field spectral directivity for Mach number profile (70) 
with a = 0.1,6 = 6, M(0) = 1 . 5 , and source position r a = 
0.75r o ,xi s = — 2.0r o for hard-wall duct (solid) and soft-wall duct 
with £ = (1, -1) (dashed), (1/2, -1) (dotted), (2, -1) (dot-dashed), 
(1,-2) (dot-dot-dashed). 

Figure - 4 Far-field spectral directivity for Mach number profile (70) 
with a = 0.1,6 = 6, M(0) = 1.5, and source position r s = 
0.75ro, £i s = — 0.5ro for hard-wall duct (solid) and soft-wall duct 
with £ = (1,-1) (dashed). 

Figure - 5 Far-field spectral directivity for Mach number profile (70) 
with a = 0.1,6 = 6, M(0) = 1.5, and source position r 3 = 
O.5r 0 ,xi s = -0.5 r 0 for hard-wall duct (solid) and soft-wall duct 
with £ = (1,-1) (dashed). 

Figure - 6 Far-field spectral directivity for Mach number profile (70) 
with a = 0.1,6 = 6, M(0) = 0.9, and source position r, = 
O.Sro,#! 3 = — 2.0r o for hard-wall duct (solid) and soft-wall duct 
with £ = (1,-1) (dashed). 



